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Abstract. Strong gravitational lensing by irregular mass distributions, such as galaxy clusters, is generally not well quantified 
by cross sections of analytic mass models. Computationally expensive ray-tracing methods have so far been necessary for 
accurate cross-section calculations. We describe a fast, semi-analytic method here which is based on surface integrals over 
high-magnification regions in the lens plane and demonstrate that it yields reliable cross sections even for complex, asymmetric 
mass distributions. The method is faster than ray-tracing simulations by factors of ~ 30 and thus suitable for large cosmological 
simulations, saving large amounts of computing time. We apply this method to a sample of galaxy cluster-sized dark matter 
haloes with simulated merger trees and show that cluster mergers approximately double the strong-lensing optical depth for 
lens redshifts zi > 0.5 and sources near Zs = 2. We believe that this result hints at one possibility for understanding the recently 
detected high arcs abundance in clusters at moderate and high redshifts, and is thus worth further studies. 



1. Introduction 

Strong gravitational lensing by galaxy clusters is a highly non- 
linear effect which is very sensitive to the details of the lens- 
ing mass distribution. The cluster core densities, the asymme- 
tries of their mass distribution, their substructures and their 
close neighbourhood all contribute to their lensing proper- 
ties. The ongoing debate about whether the observed statis- 
tics of arcs is or is not compatible with the expectations in 
the standard ACDM cosmology shows that it is not suffi- 
ciently understood yet what aspects of the source and cluster 
populations as a whole determiiies the statistics of its strong- 
lens ing effects (Bartelmann et al. 1998; Meneghetti et al. 20001 
2003b.a.: .Wambsganss et al...2004: .Dalai et aL.2004; .Li et al l 
2005tlHennawi et alJ2005r 



It is an obstacle for theoretical as well as observational 
studies that cross sections of galaxy clusters for strong lens- 
ing are costly to compute. So far, they require highly-resolved 
simulations tracing large numbers of light rays through real- 
istically simulated cluster mass distributions, used for finding 
the images of sources which need to be classified automati- 
cally (Mi ralda-Escude 1993; Bartetoann & Weiss 1994). This 
needs to be done often, i.e. for different cosmological mod- 
els and for many clusters in large cosmological volumes seen 
under many different angles, for the results to reach a reliable 
level. The recent finding that the enhanced tidal field around 
merging clusters substantially enhances strong-lensing cross 
sections (iTorri et aI-2004.) adds the necessity to study clusters 



with a time resolution which is sufficiently fine to resolve clus- 
ter merger events. 

The increasing demands to be met by reliable strong- 
lensing calculations and the wish to carry them out for vary- 
ing cosmological models calls for a substantially faster method 
than ray-tracing which should be equally reliable. We develop 
such a method here. It is based on the fact that highly elongated 
arcs occur near the critical curves in the lens plane, and that 
imaging properties near critical curves can be summarised by 
the eigenvalues of the Jacobian matrix of the lens mapping (see 
ISchneider et al. ( 1992) and Sect. 2 for detail). This allows the 
reduction of the cross-section calculation to an area integral to 
be carried out on the lens plane itself. In that sense, the method 
is analytic, but the iiTeguIar shapes of the integration domains 
require it to be carried out numerically. Since the eigenval- 
ues of the Jacobian matrix ideally describe imaging properties 
for infinitesimally small, circular sources, the method needs to 
be supplemented by techniques for taking extended, elliptical 
sources into account without losing computational speed. 

The paper describes this method and its testing. After a 
brief summary of gravitational lensing in Sect. 2, we develop 
the new method and compare it to ray-tracing techniques in 
Sect. 3. We then describe in Sect. 4 its application to samples 
of cluster-sized halos whose evolution is described by merger 
trees. Results on cluster cross sections and optical depths are 
given in Sect. 5, and Sect. 6 provides summary and discussion. 
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2. Basic Gravitational Lensing 

We briefly compile the basic equations of gravitational lens- 
ing necessary for quantifying the statistics of strong lens- 
ing events. Adopting the thin-lens approximation, the lens- 
ing mass distribution is projected perpendicular to the line of 
sight onto the lens plane, on which the physical coordinates 
^ = (^1,^2) are introduced. The coordinates on the source 
plane are rj — {rji ,172). Throughout this paper, we shall assume 
for simplicity that all sources are at the same redshift because 
their distribution in redshift is not relevant for developing our 
method. By means of the length scales ^0 and 770 = {Ds/D\)£,o 
in the lens and source planes, respectively, we can introduce 
the dimension-less coordinates x ^ and y — rj/rjo- The 
angular-diameter distances from the observer to the lens, the 
source, and from the lens to the source are Z)i s js, respectively. 

Source and image coordinates are related by the lens equa- 
tion. 



y- 



- a{x) , 



(1) 



where a{x) is the deflection angle, which is the gradient of the 
scalar lensing potential \l/{x) (see, e.g. .Schneider et al...l992) . 
The Poisson equation 



A\j/{x) = 2k{x) = 2 



E(x) 



(2) 



relates the lensing potential to the (dimension-less) conver- 
gence K, which is the surface-mass density E in units of its 
critical value Ecr. 

The lens equation Q defines a mapping from the lens to the 
source plane whose Jacobian matrix (x) describes how the 
shape of an infinitesimal source located at position y is changed 
by the lens. It can be decomposed as 



1-K-Yl -72 

-72 1-K-+71 



(3) 



where the shear components 71 2 form a trace-less tensor quan- 
tifying the lensing distortion. The magnification of an image 
relative to the source is given by the inverse Jacobian determi- 
nant. 



1 



1 



det£/{x) [1 - K-(x)]2 - ;)f (x) - 7|(jc) 



(4) 



Since the deflection angle is a gradient field, jz/ (x) is symmet- 
ric at any point on the lens plane (at least in the single-plane 
approximation which we are considering here), thus a coordi- 
nate rotation can always be found which diagonalises £/. Its 
two eigenvalues 



A, = l-K--7, Ar=l-K-+7 



(5) 



describe the distortion of an infinitesimal source in the tangen- 
tial and radial directions relative to the lens' centre and are thus 
called tangential and radial eigenvalues. At points in the lens 
plane where at least one of these two eigenvalues vanishes, the 
lens mapping becomes singular, yielding formally infinite mag- 
nifications for point sources. The set of all such points forms 
closed curves on the lens plane, the critical curves, whose im- 
ages in the source plane are the caustic curves or caustics. 



The central quantity describing the statistics of strongly dis- 
torted images with properties p is the lensing cross section Gp, 
defined as the area of the region on the source plane where 
a source with given characteristics has to lie in order to pro- 
duce at least one image with properties p. The rest of this pa- 
per deals with images having length-to-width ratios exceeding 
a fixed threshold d. The respective cross section will be denoted 
as Gd. 

This lensing cross section quantifies the probability for a 
set of sources to be lensed by a given matter distribution into 
pronounced arc-like images. If we want to know how many arcs 
are expected to form on the entire sky (that is what we observe 
or can extrapolate from observations), we must sum the cross 
sections of all individual suitable lensing masses between us 
and the source plane, and then allow for variations of the source 
redshift. The lensing optical depth for sources at redshift Zs is 



1 



4;rD2 



JQ 



N{M,z)adiM,z)dMdz, 



(6) 



where N{M,z)dz is the total number of masses M contained 
in the shell between redshifts z and z + dz- Denoting the dif- 
ferential niass function by n(M,z) jPress & Schechte3ll974^ 
Bond et al JI1 99lt ll^acev & Colelll 99/iL I1 994 ISheth & TormenI 



19991 12OO2 ') and the cosmic volume enclosed by a sphere of 
"radius" z by V{z), we can write 

dV{z) 



N{M,z) 



dz 



n{M,z) ■ 



If the number of sources with redshift between Zs and Zs 
is n^{zf)dzs, the total number of gravitational arcs is 



Nd = / ns(zs)Trf(zs)dZs 
Jo 



(7) 



dzs 



(8) 



In particular, if all sources are at the same redshift Zs, the pre- 
vious equation reduces to 



(9) 



It is important for our discussion that the lensing cross sec- 
tion a{M,z) of a single structure depends sensitively on the 
source distribution and on the internal properties of the struc- 
ture itself, in particular the slope of the inner density profile, 
the lumpiness and the asymmetries of the matter distribution 
( Bartelmann et al..,1995.. 1998: Meneghetti et al,.2003b) . These 
peculiarities of galaxy clusters are intimately Unked with the 
formation and evolution histories of the respective dark-matter 
haloes, which in turn depend on cosmology (especially on the 
matter density parameter and the cosmological constant or dark 
energy). Moreover, the spatial volume per unit redshift, the 
source populat ion and the cluster m ass function also depend 
on cosmology llLacev & Colell993l) . 

Given the importance of asymmetries in cluster mass dis- 
tributions it is not promising to search an analytic expression 
for the function (7(M,z), unless we consider high ly symmetric 
(and thus unrealistic, see Meneghetti et al. ■2003b') lens models. 
The alternative that we shall discuss in the next subsections is 
to use numerical or semi-analytic algorithms for calculating the 
lensing cross section of a set of model clusters that appropri- 
ately sample the mass and redshift range under consideration. 
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3. Lensing Cross Sections 

So far, the most reliable method for calculating cross sections 
for long and thin arcs has been using fully numerical ray- 
tracing simulations. If performed adequately, such simulations 
return realistic cross sections, but with the disadvantage of be- 
ing very expensive in terms of computational time. In view of 
cosmological applications, the mass and redshift ranges to be 
covered are large, in particular because the temporal sampling 
needs to be dense for properly resolving cluster merger events. 
The number of cross sections to be calculated can thus be 
very large. A reliable alternative to the costly ray-tracing sim- 
ulations is therefore needed. We develop here a semi-analytic 
method which fairly reproduces fully numerical lensing cross 
sections, but at a computational cost which is lowered by fac- 
tors of ^ 30. We believe this method to provide an elegant al- 
ternative to ray-tracing simulations for many applications of 
lensing statistics. 

In the following subsections, we shall first describe the 
fuUy numerical method for reference, and then the semi- 
analytic method. 



3. 1. Ray-Tracing Simulations 

Given the lensing mass distribution (which can be given as ei- 
ther an analytic density profile or a simulated density map) and 
the statistical properties of the source sample, strong-lensing 
cross sections are commonly estimated using fully numerical 
ray-tracing simulations. The method we use here was first de- 
scribed by |Miralda-Escude (Il993h and subsequently further 
devel oped and adapted to asymmetric, numerical lens mod- 
els by lBartelmann & W eissI 01994) . It has been widely used by 
iMeneghe tti et alJ l l2Q0(ll2003bl) and with several modifications 
bv IPuchwein et alJ ( l2005l) . We only address its principal fea- 
tures here, referring the interested reader to the cited papers 
and the references given therein for further detail. 



Briefly, a bundle of n x n light rays (n — 2048 here) with an 
opening angle j3 is sent from the observer through the lens. The 
opening angle depends on the lens' properties and the distances 
involved, and must be large enough to encompass the entire re- 
gion on the lens plane where strong lensing events may occur 
The deflection angle is calculated from the surface-mass distri- 
bution at all points where light rays intersect the lens plane, 
thus allowing each ray to be propagated back to the source 
plane by means of the lens equation Q. 

A set of sources is then placed on a regular and adaptive 
grid on the source plane. The sources are modelled as ellipses 
whose orientation angles and axis ratios (minor to major) are 
randomly drawn from the intervals [0,;r] and [0.5,1], respec- 
tively, with the prescription each source has the area of a circle 
with 1 arc second diameter The source-grid resolution is it- 
eratively increased near the caustics, i.e. where the magnifica- 
tion on the source plane is highest and undergoes rapid varia- 
tions. This artificial increase in the probability of strong lensing 
events, must be corrected when calculating cross sections. We 
do so by assigning a statistical weight to each source which is 
proportional to the area of the grid cell which it represents. 



The images of each source are found by identifying all rays 
of the bundle falling into the source. Simple geometrical shapes 
(ellipses, rectangles, and rings) are then fit to all images, and 
their characteristics (length, width, curvature radius, etc.) are 
determined. When an image has the property we are interested 
in (i.e. a length-to-width ratio equal to or greater than d), we 
increment the cross section by the pixel area of the source grid, 
times its statistical weight. 

3.2. Semi-Analytic Method 
3.2.1. Point Sources 

We aim at determining cross sections for the formation of 
gravitational arcs with length-to-width ratio exceeding some 
fixed threshold d. Initially assuming infinitesimal (or point- 
like) sources, the discussion in Sect. 2 implies that such arcs 
will form where the ratio 



Rix) 



At(x) 



Xt{x) 



Ar(jl 



between the eigenvalues of the lens mapping satisfies 
R{x) > d . 



(10) 



(11) 



We denote this region by Z?/ — Bi{d). By means of the lens 
equation, B; can be mapped onto an equivalent region By = 
Bs{d) on the source plane, whose area is by definition (see 
again Sect. 2) the cross section Gd we are searching for. Thus, 



(7,1= d^Tj = 77o / d^y ■ 



(12) 



The lens equation maps the infinitesimal area element on the 
lens plane to that on the source plane by means of the Jacobian 
determinant detx/ (x), thus 

a,i = rii ( |det^(x)|d2x = ri^ f , (13) 

Jb, Jb, \h[x)\ 

where Eq. ({4} was used. 

3.2.2. Extended Circular Sources 

Although this line of reasoning is exact, it fails in reproduc- 
ing simulated cross sections, for it does not account for real 
sources (and also the sources used in ray-tracing simulations) 
not being point-like. Extended sources are much more likely 
to produce strongly distorted images than point-like sources, 
and their imaging properties will only be approximated by the 
eigenvalue ratio R{x). This implies that the integral in il3\ 
has to be extended beyond the region B/, because an extended 
source can produce a long and thin arc even if it is centred 
where the eigenvalue ratio is less than d. 

Introducing extended sources into the framework just out- 
lined, we wish to convolve the eigenvalue ratio R{x) on the 
lens plane with a suitable function to be determined which sig- 
nificantly differs from zero only on the image of an extended 
source. For consistency with ray-tracing simulations, we as- 
sume circular sources with radius ^ = 0.5" for now. Effects of 
source ellipticities will be included later. 
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The obvious problem with this approach is that the prop- 
erties of the image vary wildly across the lens plane, thus the 
convolution should use a function which rapidly changes shape 
and extent across the lens plane. We can avoid this problem by 
transferring the calculations to the source plane. We first intro- 
duce the eigenvalue ratio on the source plane. Since multiple 
points Xi on the lens plane may be mapped onto a single point 
y on the source plane, we define it as R{y) = max{/?[3'(x,)]}. 
This ratio R{y) is then convolved with a function g{y) which 
differs significantly from zero only on the circular area covered 
by a source, which is the same everywhere on the source plane. 
Let the convolution in the source plane be h{y) = R{y) *g{y), 
then we put h{x) — h[y{x)] to obtain the convolved function on 
the lens plane, as desired. Again, we need to take into account 
that single points on the source plane may be mapped on multi- 
ple points in the lens plane. Assigning R{x) to R[y{x)] now sets 
the convolved length-to-width ratio equal on all image points, 
which is not exact, but a tolerable error The problem is now re- 
duced to carrying out the convolution h{y), which is discussed 
in the Appendix. As described there, we speed up the convo- 
lution by approximating it by a simple multiplication after es- 
sentially assuming that the eigenvalues of the lens mapping and 
their ratio do not change significantly across a single source. 

Finally, what form should be chosen for the function ^(y)? 
A two-dimensional Gaussian of width ^ may appear intuitive, 
but the ray-tracing simulations we use for reference do not 
adopt a surface-brightness profile for sources because they sim- 
ply bundle all rays falling into the ellipse representing a source. 
Thus, a choice for g consistent with the ray-tracing simulations 
is a step function with width ^, 



else 



where y^3§y < 1 



(14) 



where ^ is a matrix describing the shape of the source. Since 
we are only considering circular sources of radius ^, ^ has 
the form ^ = /C^, where is the unit matrix. The factor 
(tti^^)^' in J14t normalises g{y) to unity. Substituting li{x) for 
R{x), we can again apply Eq. il 1> to obtain a new region B/ 
that will now of course be larger than before. 

3.2.3. Extended Elliptical Sources 

We finally have to account for elliptical rather than circu- 
lar source s. It is quite obvious and was shown by many au- 
thors (e.g. lBartelmann et alJll995UKeetonl2001 *) that elliptical 
sources are more likely to produce strongly distorted images. 
We adopt the simple and elegant formalism by iKeetorJ 



( 1200 11) . He showed that an elliptical source at position y{x) with 
axis ratio qs=a/b — qs\y{x)] and orientation angle Q — Q[y{x)] 
is imaged as an alUpse with an axis ratio of 



9obs = 



lT + \/T^--AD 



(15) 



where T and D are the trace and the determinant of the matrix 
that describes the image ellipse. They can be expressed as func- 



tions of the intrinsic source properties and of the (convolved) 
lensing distortion: 

T = h- + q^. + {h^^l){q^.- I) cos^e , D = h^q]. (16) 

Again, h{x) can be replaced by ^obs to obtain a modification 
of the region B/, again larger than before, whose area is now a 
close approximation to the cross section we seek to determine, 
accounting for extended and elliptical sources. 

3.2.4. Comparison to Ray-Tracing 

As a final step, we test the accuracy of our calculations and ap- 
proximations, especially those concerning the replacement of 
the convolution by a multiplication (cf. the Appendix), by com- 
paring lensing cross sections obtained with the semi-analytic 
method to their fully numerical counterparts obtained by ray- 
tracing. Since the main purpose of this work is to estimate the 
effects of cluster mergers, we compare the cross sections as 
they evolve while two dark matte r haloes merge. The halos 
are rnodelled as NFW profi les (see iNavarro et alJI 19951 Il996l 
ll997t F ole & Lacev * 1996) whose lensing potential is ellipti- 
cally distorted to have ellipticity e = 0.3 (see Sect. 4 for detail). 
The results of the comparison are shown in Fig. [l] for various 
masses, redshifts and thresholds d. 

Figure[OwaiTants one observation. The ray-tracing simula- 
tion was repeated five times for each modelled merger event, 
each time changing the seed for the random generation of 
source ellipticities and orientations. As the figure shows, this 
has quite a significant effect on the numerical cross sections, 
causing a substantial scatter. The reason can be easily under- 
stood. When a given source produces an arc, variations in the 
source's intrinsic ellipticity and in the orientation with respect 
to the caustic structure can easily push the length-to-width ratio 
of the image arc below or above the chosen threshold d. This 
causes the irregular trend seen in Figure ^ It is worth noting 
that the same problem in principle also affects the semi-analytic 
results. In this case, however, ellipticity and orientation are as- 
signed to each ray traced back through the lens plane simply to 
identify the the extended region over which we integrate. Since 
the number of rays is much larger than the number of individ- 
ual sources used in the fully numerical algorithm, the random 
scatter in the semi-analytic results is very small. In fact, the 
fluctuations are of the order of the width of the black curve and 
thus omitted. 

Moreover, it must be noted that the results given by the ray- 
tracing code we use might differ slightly from other codes us- 
ing different resolution, different image finding algorithms and 
different ways to fit the image shapes. 

Analysing Fig.^ we can draw two interesting conclusions. 
The first is that, reassuringly, the numerical and the semi- 
analytic cross section agree excellently. This means that the ap- 
proximations we made are substantially correct. There is, how- 
ever, a small discrepancy for higher-redshift sources, shown in 
the lower panels. We believe that this discrepancy is due to 
the fact that the scale over which the lensing properties on the 
source plane change significantly is in this case comparable to 
the source dimension, thus one of the approximations we made 
in the Appendix fails, namely that the function /^[^(.x)] does 
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Fig. 1. Lensing cross sections for long and thin arcs with length-to-width ratio exceeding c/ as a function of the distance between 
the centres of two dark matter haloes, for two (lens and source) redshifts. The haloes are modelled as NFW density profiles with 
elliptically distorted lensing potential with ellipticity e = 0.3 (see Sect. 4). Green lines show the results of ray-tracing simulations 
for five different realisations of the source distribution. Black lines are the results obtained with our semi-analytic method. 



not vary much across a source. Nonetheless, the discrepancy 
is nowhere larger than 20%, which is more than acceptable for 
our purposes, in particular in view of the considerable scatter 
in the ray- tracing results. 

The second observation is that the behaviour of the cross 
sections (both numerical and semi-analytic) as a function of 
the distance betw een the haloes closely reflects that found by 
Irbrri et aTl ( 120041) . These authors studied the effect of a single 
merger process on the strong-lensing efficiency of a numeri- 
cally simulated dark matter halo, using ray-tracing as described 
in Sect. 3. They found that, while the substructure is swallowed 
by the main halo, there is a first peak in the lensing cross sec- 
tion when the increasing shear field between the lumps causes 
the critical lines to merge, a local minimum while cusps dis- 
appear in the caustic branches, and another peak caused by the 
increased convergence when the two density profiles overlap. 
All these features are recovered in the panels of the Fig.^ 

This level of agreement between ray-tracing and semi- 
analytic cross sections show that the semi-analytic method de- 
veloped above for cross-section calculations is essentially cor- 
rect and constitutes a valid and useful, ~ 30 times faster alter- 
native to the costly ray-tracing simulations. 



4. Lensing Optical Depthi 

Given reliable lensing cross sections, lensing optical depths 
need to be determined following Eq. (|6j. For that purpose, we 
use the merger trees for a set of 46 numerically simulated dark- 
matter haloes, whose main properties are described in the next 
subsection. The merger tree of a dark-matter halo provides two 
types of information. First, we know how the mass of each halo 
evolves with time or redshift. Second, we know at which red- 
shift merger events happen with substructures of known mass. 
Regarding this, it is worth recalling that the evolution of a dark- 
matter halo is characterised by the continuous accretion of in- 
falling external material. We account for afl mergers in which 
the main halo accretes a sub-haloes with at least 5% of the main 
halo's mass. Nonetheless, it should be kept in mind that halos 
continuously accrete matter apart from merger events. 

4.1. Halo Model 

We base our study on the merger trees on a sample of 46 nu- 
merically simulated dark-matter haloes. These haloes were re- 
simulations at higher resolution of a large-scale cosmological 
simulation. The cosmological model used was a "concordance" 
ACDM model, having a cosmological constant of CIa.o = 0.7, 
a (dark) matter density of ii^.o =0.3 and a Hubble constant 
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Table 1. Present-day masses and virial radii for the three most 
massive and the three least massive haloes in our sample. 



Halo Id. 


Virial Mass 


Virial Radius 






[/j-'Mpc] 


g8-a 


2.289 10'^ 


2.146 


gl-a 


1.530 10'^ 


1.876 


g72-a 


1.374 10'^ 


1.810 


g696-y 


5.219 10^-* 


0.608 


g696-z 


5.171 10^^ 


0.607 


g696-# 


5.060 10'^ 


0.602 



of h ~ 0.7. The power spectrum of the primordial density- 
fluctuation field is scale invariant (i.e. with a spectral index of 
n ~ 1), and the nns linear density fluctuations on a comoving 
scale of Mpc is ag = 0.9, which is the typical value re- 
quired to match the local abundance of massive galaxy clusters 
(White et al. 1993; Eke et al. 1996). The mass of dark-matter 
particles is w = 1.3 x lO^/i^^M©. In Tab.[2we list the present 
masses and virial radii for the three most massive and the tree 
least massive haloes in the sample. 

The present masses of the halo models vary between about 
5 X IO'^/i^'Mq (just exceeding the mass of a massive cD 
galaxy) and more than 2 x 10^^ h^^ Mq, which is characteristic 
of a rich galaxy cluster. 

We note in passing that, owing to the dependence of the 
cluster-evolution history on cosmology, the typical redshift of 
structure formation in the cosmological model used here is 
higher than in an Einstein-de Sitter universe, and lower than in 
an open low-density universe. Models of dark energy alterna- 
tive to the cosmological constant typically shift this redshift to- 
wards higher redshifts, thus changing the contribution of clus- 
ter mergers to arc statistics. We are planning a future study to 
address this interesting property. 

For describing specifically each individual merger process 
and its effects on gravitational-arc statistics, we model each 
dark-matter halo as NFW spheres. Their density profile is 



r/r,(l+r/rs)2 ' 



(17) 



which fits the density profiles of numerically simulated dark- 
matter haloes and fairly describes dark structures over a wide 
range of masses and in different cosmologies. The scale density 
Ps and the scale radius are linked by the halo concentration, 
c = '"200 /'"si where r2oo is the radius enclosing a mean density 
of 200 times the closure density of the Universe. As a gen- 
eral trend, the concentration reflects the mean density of the 
Universe at the collapse redshift of the halo, thus it is higher 
for lower masses. Several ways to quantitatively relate the con- 
centration to the halo mass have been proposed. We follow the 
prescription of Eke et al. ( 2001 )- The analytic fit In} allows 
direct calculations of the main lensing properties. In particular, 
the projected density profile for a NFW lens is 



and its deflection angle is 

a(x)=4psrs^, 

X 

where the functions f{x) and g{x) are 



(19) 




(20) 



nx) 



2psrs^ : 



(18) 



dBartelmanrJ 1 1 99^) . Here and in the remainder of this sub- 
section, X is the dimensionless distance from the lens' cen- 
tre. We additionally adopt elliptically distorted lensing poten- 
tials. Elliptical potentials are not equivalent to elliptical den- 
sity distributions, but allow for much more straightforward cal- 
culations (Schnsideret alg, 1992). and the results are equally 
vahd fo r our purposes (joolse & Kneibl2002HMeneghetti et alJ 
!2003b). In their work.'M eneghetti et alJ (^2003b') also estimated 
the value of the ellipticity e for the lensing potential that pro- 
duces the best fit to the deflection-angle maps of simulated 
haloes. Following their result, we adopt e = 0.3 throughout. 



4.2. Modelling Halo Mergers 

Once the redshift is fixed, the halo we consider may be iso- 
lated or in interaction with a substructure. In the first case the 
deflection angle can be obtained directly by means of equation 
( I19> . In the second, we can sum the deflection angles of the 
two structures at each point on the lens Pl™s, owing to the lin - 
earity of the problem (Schneider et al. 1992; Torri et al. '2004V 
In both cases, the eigenvalues of the local mapping follow from 
the deflection-angle maps by differentiation, and can be applied 
in the semi-analytic method for calculating cross sections as 
described in Sect. 3. 

Individual merger events are modelled as follows: The main 
halo and the substructure start at a mutual distance equal to the 
sum of their virial radii. Then the centres approach at a con- 
stant speed, calculated by the ratio of the initial dist ance to the 
typical time scale of a merger event, i.e. 0.9 Gyr ("T orri et aP 
12004; Tormen et al. 2004) . The process is assumed to be con- 
cluded when the two density profiles overlap perfectly. The di- 
rection of approach, always assumed to be perpendicular to the 
observer's line of sight, is parallel to the directions of the ma- 
jor axes of the lensing potential, which are also assumed to be 
paralle l. This last assumption is well justified by recent work 
dUee et al. 2005: Hopkins et aL 2005.) pointing out, both with 
analytic models and numerical analyses, that there is an intrin- 
sic alignment between a dark-matter halo and the suiTound- 
ing haloes and sub-lumps due to the tidal field of the major 
halo itself. We note that assuming all mergers to proceed along 
directions perpendicular to the line-of-sight slightly underes- 
timates the cumulative contribution of cluster mergers to the 
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lensing optical depth, because if the merger direction is par- 
tially aligned with the line-of-sight, the time spent by the sys- 
tem in configurations leading to high cross sections is longer 
(see again Fig. |3jl. Several tests we made show however that 
this underestimation is typically of order 20%, thus we neglect 
it here. 

Once we know for each halo at every redshift the cross sec- 
tion with and without the effect of mergers, we can compute 
lensing optical depths. Since we sample the mass range dis- 
cretely, we cannot apply Eq. (|6} directly, but need to approxi- 
mate it as 



n-1 
1=1 



Mi 



N{M,z)dM 



dz 



(21) 



where the masses Mj (1 < i ^ n) have to be sorted from the 
smallest to the largest at each redshift step, and the quantity 
dd.i{z) is defined as 



OdAz) ^ \ [OdiMi.z) + ad{Mi+uz)] 



(22) 



This is essentially equivalent to assigning to all structures 
with mass within [M,-,M(+i] the average cross section of the 
haloes with masses M, and M,+i, weighted with the number 
density of massive structures within that interval. 



46 halos which still contributes to strong lensing is shown in 
Fig. 12] as a function of lens redshift with fixed source redshift 
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4.3. Mass Cut-Off 

Since a sufficient condition for strong lensing is satisfied if the 
surface density of a lens exceeds the critical density some- 
where, any lens model with a cuspy density profile such as 
the NFW profile will formally be a strong lens and thus pro- 
duce critical curves and caustics. However, the caustics of low- 
mass halos will be smaller than the typical source galaxies. 
Averaging the local distortion due to a low-mass lens across 
an extended source will then lead to a small or negligible to- 
tal distortion. This implies that the mass necessary for halos to 
cause large arcs is bounded from below by the requirement that 
the halo's caustic must be sufficiently larger than the available 
sources. 

This mass limit obviously depends on the lens and source 
redshifts due to the geometrical sensitivity of lensing. Higher 
masses are required at low and high redshifts for lensing effects 
comparable to lenses at intermediate redshifts. In addition, the 
caustic structures can change substantially during major halo 
mergers. As a sub-halo approaches the main halo, the initially 
separated caustics of the two merging components will increase 
and merge to form a larger caustic. Thus, even though the to- 
tal mass is unchanged, the mass limit for strong lensing may 
decrease while a merger proceeds. Even halos which are indi- 
vidually not massive enough for arc formation may be pushed 
above the mass limit while they merge. In view of the exponen- 
tially dropping cluster mass function, this is potentially a huge 
effect. 

Thus, we have to monitor the extent of the caustics as we 
compute the lensing optical depth of a halo sample, taking into 
account that the mass limit may change rapidly as halos merge 
with sub-halos. The lowest mass of a halo from our sample of 



Fig. 2. Mass of the lowest-mass halo producing large arcs in the 
sample of 46 halos which we use here to compute the lensing 
optical depth. The source redshift is Zs = 2. The overall trend in 
the curve reflects the geometrical lensing sensitivity, while the 
fluctuations and the depth of the minimum reflect that merg- 
ers can lift low-mass halos above the strong-lensing threshold 
which would otherwise not be capable of strong lensing. 

Following these prescriptions, we are able to calculate the 
strong-lensing cross section of each model halo at every red- 
shift step of the simulation, first ignoring the effects of merger 
processes, then accounting for them. Thus, we sort the masses 
from the smallest to the largest at every redshift step and cal- 
culate the quantities (Jd.iiz)- Finally, we calculate the lensing 
optical depth with and without the effect of merger processes. 

5. Results 

5. 1. Cross Sections 

We calculate the optical depth for the formation of gravitational 
arcs with length-to-width ratio grater than an arbitrary thresh- 
old d. We used two choices, d — 7.5 and d — 10. 

Before discussing the calculation of the optical depth it is 
interesting and useful to study the behaviour of lensing cross 
sections of individual haloes with redshift. In Fig.|3]we show 
the cross sections for four of the most massive haloes, both with 
d — 7.5 (top panels) and d — 10 (bottom panels). The sources 
are put at a fixed redshift z, = 2. Apart from the obvious fact 
that cross sections for arcs with higher length-to-width ratio 
are smaller than those including shorter arcs, we see that all 
cross sections tend to zero when the redshift approaches zero 
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Fig. 3. Evolution of the lensing cross section for gravitational arcs with length-to-width ratio equal to or greater than d = 1 .5 (top 
panels) and c/ = 10 (bottom panels) for four of the most massive haloes in the sample. Sources are at redshift z% = 2. Red-dashed 
lines show cross sections calculated without taking account of merger processes. Black lines show cross sections enhanced by 
cluster interactions. Filled blue dots are the counterparts of the black lines obtained from fully numerical ray-tracing simulations. 



or the source redshift. This "suppression" due to the lensing ef- 
ficiency is caused by the geometry of the problem, and in par- 
ticular by the fact that lenses very close to the sources or the ob- 
server have arbitrarily large critical density. Moreover we note 
that the increase in the cross section due to merger events can 
exceed half an order of magnitude in some cases. Torri et al. 
(120041) found an increase in the lensing cross section up to an 
order of magnitude, but here, even if mergers happen at red- 
shifts with higher lensing efficiency {z\ ~ 0.2 — 0.5), the masses 
involved are lower, so we cannot reach such higher increases. In 
particular, high-redshift mergers are quite inefficient in boost- 
ing total lensing cross sections, both because of low involved 
masses and proximity to the source plane. To further test the 
reliability of our semi-analytic calculations we show as filled 
blue dots in Fig.|3]the cross sections obtained from ray-tracing 
simulations. The agreement is again reassuringly good. 

5.2. Optical Depths 

In Fig. 13 we show the optical depth per unit redshift, i.e. the 
contribution to the lensing optical depth given by structures in 
different redshift bins per unit lens redshift. In other words, 
the integrand of the redshift integral in Eq. i21\ is shown. The 
integration over mass is carried out above the mass limit illus- 
trated in Fig. 12] i.e. all halos are included which have caus- 
tic structures sufficiently larger than individual source galax- 
ies. The two panels are for d = 7.5 and d ~ 10, respectively. 
Four curves are shown in each panel, obtained ignoring merg- 
ers (smooth, dashed curves) and taking mergers into account 
(solid curves). In both panels, the upper and lower curves refer 
to source redshifts Zs = 2 and Zs = 1, respectively. 



The overall trend of the differential optical depth in Fig. |4] 
resembles individual cross sections, i.e. it drops to zero as the 
lenses approach the observer or the sources. The dashed (up- 
per) curve for Zs = 2 broadly peaks at redshift Zp ~ 0.4, slightly 
larger than the typical redshift for the peak of the individual 
cross sections shown in Fig. |3l This is simply due to the fact 
that in the differential optical depth the lensing cross sections 
are weighted with the number density of structures within mass 
bins. It is interesting to note that the same peak occurs even in 
the corresponding solid curve, thus it is not due to dynamical 
processes in the cluster lenses, but rather to the combination 
of the mass evolution of the lenses with the particularly high 
lensing efficiency for clusters at that redshift. 

Apart from that, the most remarkable result shown by the 
solid curves is that the impact of cluster mergers is important 
in particular at moderate and high redshifts, 0.5 < z < 0.8. The 
pronounced peaks in the differential optical depth seen there 
even after averaging over the halo sample indicate that cluster 
mergers can substantially increase the lensing optical depth of 
high-redshift clusters. Above redshift 0.5, mergers almost dou- 
ble the optical depth. 

Shifting the source plane from Zs = 2 to Zs = 1 significantly 
lowers the total optical depth as well as the impact of cluster 
mergers on the optical depth per unit redshift. The first effect 
is obviously due to the fact that for lower source redshift, the 
redshift interval of high lensing efficiency narrows. The second 
effect reflects that sources at lower redshift miss a significant 
part of the lensing haloes' formation history and merger pro- 
cesses related with it. 
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Fig. 4. The evolution with redshift of the optical depth per unit 
redshift for gravitational arcs with length-to-width ratios equal 
to or larger than d. Solid lines include the effect of cluster merg- 
ers, while dashed lines do not. Upper curves refer to sources 
with Zs—2 and lower curves to sources with Zs = 1 • 



5.3. Sources Properties 

A full analysis of the effect of various source properties on the 
cross sections of individual haloes is well beyond the scope of 
this paper Moreover, the method we have outlined in the previ- 
ous subsections is probably not the ideal tool for that purpose. 
For example, it would be of interest to check the variation of 
cross sections with the source size. Fig.|5]shows that numerical 
and semi-analytic cross sections agrees (as already shown) and 
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Fig. 5. Cross section for gravitational arcs with length-to- width 
ratio exceeding d=10 for a dark-matter halo of lO'^M©/;^' as 
a function of the source size in arc sec. Black and green lines 
show the results of the semi-analytic and of ray-tracing simula- 
tions, respectively, which are repeated with different random- 
number seeds for the ellipticities and the position angles of the 
sources, as in Fig.fD The source redshift is Zs = 1, and the lens 
redshift is zi =0.3. 



initially increase with increasing source size. However, when 
the sources become too large compared to the size of the caus- 
tic structure, the assumption that the lensing properties are ap- 
proximatively constant across the area of a source fails (see 
the Appendix for details). Thus, while numerical cross sections 
start to decrease because the sources are too large to be effi- 
ciently distorted, the semi-analytic cross section increases dra- 
matically. We note that, as in the rest of this paper, dark matter 
haloes are modelled with NFW density profiles and elliptically 
distorted lensing potential. 

Nonetheless, some testable effects can be briefly addressed 
here. The first is the influence of the shape (circular or ellipti- 
cal) of the sources on the lensing efficiency. On this ground we 
can compare our results with Keeto n ( 2001). 

Fig-Elshows the cross section of several dark-matter haloes 
of increasing mass as a function of the source ellipticity, as 
explained in the caption. We see in the Figure that the cross 
sections for small elliptical sources exceed those of small cir- 
cular sources by a factor of « 2. This agrees with the findings 
of Keeto n (2P01J- It is quite interesting to note that the corre- 
sponding increase is somewhat lower for larger sources. This 
is due to the fact that the cross section is on the average higher, 
thus the contribution from the source ellipticity is relatively less 
important. 

Another quite important effect is the change of the cross 
section with the source redshift. Investigating this, we will keep 
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Fig. 6. Cross section for gravitational arcs with length-to-width 
ratio exceeding c/ = 10 for several haloes of increasing mass. 
Black and red lines refer to point-like or extended sources with 
area equal to that of a circle of radius 0.5", respectively. Short 
dashed lines indicate circular sources, long dashed lines refer to 
sources with random ellipticity drawn from the interval [0,0.5] 
(^.s G [0.5, 1]), and solid lines refer to sources with ellipticity 
0.5. Sources are at redshift Zs = 1 and the lens redshift is z\ — 
0.3. 

the source size fixed at a radius of 0.5", since the angular di- 
ameter distance changes only a little over redshift 1 . 

In Fig.0we plot the cross section for a single dark-matter 
halo of IO^^Mq/i-i as a function of the source redshift. We 
adopt three different lens redshifts, 0.3, 0.5 and 0.7 respec- 
tively. The low-redshift end of that plot shows that the closer 
the lens is to the source, the lower the cross section is due to 
the geometrical suppression of the halo's lensing efficiency (see 
the discussion in Sect. 5). The general trend is a rapid increase 
of the lensing cross section with increasing lens redshift fol- 
lowed by an almost constant phase and a slightly decreasing 
phase. This trend agrees with the general evolution of the lens- 
ing efficiency with the source redshift. 

These few examples show how the influence of the source 
redshift and intrinsic properties on a cluster's lensing efficiency 
is significant, in agreement with earlier work, and a deeper 
analysis of the consequences is certainly needed. 

6. Summary and Discussion 

We have described a novel method for semi-analytically cal- 
culating strong-lensing cross sections of galaxy clusters. The 
method first approximates the length-to-width ratio of images 
by the ratio of the eigenvalues of the Jacobian matrix of the 
lens mapping. The requirement that this ratio exceed a fixed 



threshold defines a stripe on both sides of the (tangential) caus- 
tic curve whose area approximates the cluster's strong-lensing 
cross section. 

This approach would be valid for infinitesimally small, cir- 
cular sources. Extending it for elliptical sourc es is stra i ghtfor - 
ward using the elegant technique developed bv iKeetoiil ( l200 ll) . 
Extended sources can be taken into account after convolving 
the eigenvalue ratio with a suitable window function quantify- 
ing the source size. In order to speed up this convolution, we 
approximate it by a simple multiplication. 

We have tested this method by detailed comparison with 
cross sections obtained from full ray-tracing simulations. We 
found excellent agreement within the (considerable) error bars 
of the ray-tracing results for a variety of lens masses and lens 
and source redshifts. Deviations occur for very weak lenses 
whose caustics are so small that the crucial assumption of 
the eigenvalue ratio's not changing much across sources is no 
longer satisfied. In particular, our tests revealed that cross sec- 
tions rapidly changing during merger events are well repro- 
duced by our new method. 

We then proceeded to apply the technique to a sample of 
halos whose history is described by simulated merger trees. 
The halos themselves are modelled as pseudo-elliptical halos 
with NEW density profile whose mass is given as a function 
of redshift by the merger tree. The merger trees are obtained 
from a sample of 46 cluster-sized halos numerically simulated 
in a cosmological volume. We followed the evolution of the ha- 
los by simulating merger events at times when the merger trees 
signal the accretion of a sub-halo with mass comparable to that 
of the main halo. 
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Fig. 7. Cross section for arcs with length-to- width ratio exceed- 
ing d — 10 for a dark matter halo of IO'^Mq/z^^ as a function 
of the source redshift. Three different lens redshifts are consid- 
ered, as labeled in the plot. 
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This technique allowed us to study the total optical depth 
for arc formation by the simulated cluster sample at a time 
resolution high enough for properly following merger events. 
Comparing the results to those obtained ignoring mergers, we 
found that the arc optical depth produced by clusters with mod- 
erated and high redshifts, z > 0.5, is almost doubled by merg- 
ers. 

The resuts just outlined may be potentially relevant 
in view of the high frenquency of arcs recentl y detected 
in clusters at moderate an d high redshifts (Gladders^g^^ 
l2003t | Zaritskv & Gonzale32 003). For example, Gladders et alJ 
ll2003 ) argue hat some physical process must boost the lensing 
efficiency of high-redshift clusters which is likely connected 
with the dynamics and formation histories of the lensing clus- 
ters, asjnergerprocesses are. This conclusion was supported 
bv dHoresh et aljr2005 ) after this paper was first submitted. 
lHoreshetal.l(l2005h compare the lensing efficiency of matched 
observed and simulated (low-redshift) galaxy clusters with a 
realistic source population taken from the Hubble Deep Field. 
They find that real clusters are a Uttle bit more efficient lenses 
than simulated ones, and though the difference is marginally 
significant, they argue it could be due to a selection effect con- 
nected with cluster mergers. In fact, cluster mergers increase 
not only the lensing efficiency, but also the X-r ay luminosity 
Randall et al.l2002i) . and the real clusters used bv lHoresh et al.l 



1 20051) are X-ray selected. In other words, the observed lensing 
clusters could be a biased sub-set of the entire cluster popula- 
tion. 

The method described and developed here reduces compu- 
tation times for strong-lensing cross sections by factors of ^ 30 
compared to ray-tracing simulations. It thus becomes feasible 
to reliably compute strong-lensing probabilities describing an 
evolving cluster population by halos accreting mass as encoded 
by simulated merger trees, whose merging events can be stud- 
ied at high time resolution. We shall apply this technique in 
forthcoming studies with special emphasis on the importance 
of cluster evolution for strong lensing. 
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Appendix A: Approximate Convolution of a 
Function With a Step Function on the Source 
Plane 

Consider an arbitrary function R{y) defined on the source 
plane. Suppose we wish to convolve R{y) with another func- 
tion g{y), defined as in Eq. J14> . The convolution is 



h(y)^{R*g){y)= I R(z)g{y^ 



(A.l) 
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Without loss of generality, any given pointy on the source plane 
can be chosen as the coordinate origin. This means that y = Q, 
hence 



/2(0) = (^*^)(0)= I R{z)g{z)A^ 



(A.2) 



Now, we choose a position u on the lens plane such that z = 
z{u). Applying the lens mapping to the convolution above we 
obtain 



h{0) 



R{u)g{u) 



Im(«) 



(A.3) 



We assume that R{y) does not vary much across a source. This 
assumption is satisfied in almost all interesting cases, except 
when the sources are at high redshift and the lens is close to 
them. In that case the critical curves are very small, thus the 
typical scale on which the lensing properties vary may be com- 
parable with the angular extent of a source. However, even in 
that case the results of our method remain good, in particular 
in view of the substantial scatter in the ray-tracing results. 

Within this assumption, we can expand the function R{x) 
into a Taylor series around zero, obtaining (summing over re- 
peated indices) 



h{Q) « R{Q) / g{u) 
1 d^R{0) 



\n{u)\ 

UiUjg(u) 



dR{0) f , , d^M 
/ Uig(u) 



\^iu)\ 



(A.4) 



The first integral is unity by normalisation, and the second van- 
ishes because of the symmetry of g{x). Thus, Eq. (IA.4t reduces 
to 

h{y) « R{0) + -^--±-L i^UiUjg{u 



2 dujduj Jr 
, , 1 d^R{0) 

= /?(o) + -^^a- 

2 auiOUj 



(A.5) 



Thus, we have to carry out the three integrals iln, Q.22 and 
D.12 = ^21- We show the explicit calculation only for the first, 
as the others are quite similar: 



1 



d^u 



\^i{u)\ kC^Jd |m(« 



(A.6) 



where D is the set of all positions x on the lens plane where 
x^Tx < 1. The matrix F defines the shape of the image 
formed from the source and can be written as F = st^"^ SSs^ = 
s/'^ / 11,^ . Obviously, the eigenvalues of T are '^j^'^ and 
Aj^/^^. We can now rotate into a reference frame in which 
and thus also F are diagonal. This is achieved by a rotating 
about an angle 



1 , (Ji 
(p = - arctan — 

2 Vri 



(A.7) 



where 71 and 72 are the two shear components at the origin. 
This rotation is described by the orthogonal matrix 



cos 9 — sm^ 
sin (p cos (p 



(A.8) 



With V = Ru, this rotation, transforms Eq. iA.6i into 



£lii = [ (v?cos^9 + V2sin^(p) 

Jd 

d^v 

(2viV2Cos^sin9)—- --- 

D ImMI 

-, d\ 



dh' 



cos^ (p 



sin^ (p 



Jd MmMI ?rC2 



d^v 



(A.9) 



In the new coordinate system, D is the set of all positions v of 
the lens plane where 



(A. 10) 



Now we can introduce polar elliptical coordinates (p, 0) by 



c c 

^1 = TTif"^"^^ ' ^2 -prrP sin , 
|A,| |Ai-| 

in terms of which Eq. ( IA.9l l becomes 



(A.ll) 



i2, 



^^cos^^ 

7t Jo 

C^sin^(p 

TT Jo 



2n rl 



3 2n ^P^^ 

^ A,2(p,0) 



3 . 2q dpde 
p sm 0- 



(A.12) 



A2(p,0) 

As a second approximation, we shall assume that the tan- 
gential and radial eigenvalues do not vary much across a single 
image, so we can replace the eigenvalues by their mean values 
across the image. Then, 

C^cos^^ f^'^ '■^ 
Jo 



7r(A2) Jo 
4 



p^cos2 0dpd0 



p^sin2 0dpd0 



COs2 (p 



sin2 (p 
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Similarly, we find that 



^22 = ^ 



^2 / cos2 (p 



sin2 (p 



and 
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Substituting into Eq. (IA.5> . we obtain 

1 d^R{0) f^^cos^-(p C^sin2^ 



(A. 13) 
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+ 2^xf[~Mx^^^wr 

d^R{0) f coscpsincp ^2cos^sin^\ 
+ d^2 [ 4(A2) i(A2^ J ■ ^^-^^^ 

Within the framework of our approximations, we can thus re- 
place the value of the function R at a point of the lens plane 
with its convolution on the source plane at the corresponding 
point, represented by Eq. ( IA.16t . In this way, we account for 
finite source sizes. 



